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The polarization effect of the 31 reflections within the mirror train of an auxiliary telescope of 
the Very Large Telescope Interferometer is calculated as a function of pointing direction and rotator 
angle. With a rough estimate of the mean complex index of refraction of the reflecting surfaces, 
their Jones matrices are concatenated while tracing a ray from Ml up to the feeding mirror in front 
of an instrument. The net effect is summarized in terms of the axis ratio of the polarization ellipse 
of star light that was circularly polarized above the primary mirror. 
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I. COUDE TRAIN OF AUXILIARY 
TELESCOPES 

This manuscript is concerned with the polarimetric 
characterization of the mirrors relevant to observations 
with the Very Large Telescope Interferometer (VLTI). 
The competition for observing time with the 8-m Unit 
Telescopes implies that the 1.8-m Auxiliary Telescopes 
(AT's) are more important to interferometry from a sta- 
tistical point of view, and only these will be considered 
here. 

The results demonstrate the polarization inherent to 
one beam delivered by one telescope. The more chal- 
lenging application to compute the differential polarime- 
try differentiating two beams by two different telescopes 
[lOl is not addressed; in the sense that coherency re- 
mains intact supposing the phases in the two polarization 
sub-channels are modified by the same amount in both 
beams, interferometry is not threatened. However, if 
an interferometric instrument depends on a scheme that 
measures the fluxes in the two polarizations separately, 
imbalanced fluxes induce unfavorable signal-to-noise ra- 
tios [iil,!!^. 

The configuration of the reflecting mirror optics is out- 
lined in Figure [T] [T^ , which traces five rays at various 
distances to the optical axis starting with the first reflec- 
tion off the 1.8-m primary Ml: 

• 8 reflections in the Coude train. (The design of 
the unit telescopes ^ up to M8 is similar. The 
first three of these appear to be understood [46j.) 
Incidence angles at M3 and M4 are typically close 
to 45°; 

• 3 reflections (Dl, D2, D3) inside a star rotator fl3| 
with incidence angles of 19° and two times 54.5°; 

• and 8 reflections M9 to M16 inside the star separa- 
tor [l2|, incidence at M9 close to 45° and otherwise 
< 16°. 
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After these first 19 reflections, the light beam leaves M16 
and heads with a nominal diameter of 80 mm into the 
light duct. The followup reflections are included tak- 
ing rough information on the geometries (conic constants 
and apex positions) with a ruler from a variety of design 
sketches [li,[23,|4|: 

• 1 reflection at 45° off what was M12 in the nomen- 
clature before introduction of the star separator 
PH to align the beam with the ±J7 axis of the 
delay line tunnel; 

• 5 reflections inside the main delay line cat's eye [23| : 
incidence typically < 5°; 

• 1 reflection at 45° pushing the beam from the delay 
line tunnel into the VLTI laboratory, into the +V 
direction, with the aid of M16 in the old nomencla- 
ture [3^; 

• 3 reflections below 4° inside the beam compressor, 
quenching the beam diameter to 18 mm; 

• 1 reflection in the "switch yard" with incidence an- 
gle 45° either forward to the differential delay lines 
or to the instruments (The 5 reffections inside dif- 
ferential delay lines are only relevant to a mode of 
one instrument and not included for that reason 

SEl); 

• 1 dichroic reflection at 45° into the -\-V direction; 
feeding individual instruments. 

Numbers in Section IIIII are a snapshot including these 
31 reflections, deflning the interface towards the first op- 
tical element of a generic VLTI instrument. 

Figure [3] clarifies the azimuth definitions of this 
manuscript in a projection looking from the zenith above 
the telescope down the vertical axis onto M3 or M8. The 
position of the star fixes the angle A and zenith angle z 
in local horizontal V , U, W coordinates. The rotator an- 
gle r (which rotates the image by 2r) is defined with the 
same azimuth reference, characterized by the position of 
D2. 
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FIG. 1: The AT mirror cascade. Configurations with two 
different iocations of the Nasmyth platform, which place mir- 
rors M4 to M7 at opposite sides of the vertical azimuth axis, 
are shown in red (N-|-) and in green (N-). The magnified 
view below is painted with one color only; in this section the 
geometric paths of N-l- and N- do not differ. 




FIG. 2: The coordinates in a bird's view, looking down the 
negative W direction of Figure[T] M8, Dl, D3 and M9 appear 
stacked below M3. The azimuth A is defined relative to the 
direction in which the beam leaves the Coude train into the 
duct. The two prospective directions from M3 to the Nasmyth 
focus and to M4 co-rotate with the vertical axis, A + 90° for 
N+, A - 90° for N-. 



II. RAY TRACING 
A. Reflective Surface Ellipsometry 

Given the incidence angles on the mirror surfaces, the 
polarimetric effects on the beam can be calculated "in 
principle" (to quote Tinbergen [31 ) if the complex valued 
refractive index n of these are available at the wavelength 
of interest. 

The complex- valued index of refraction for Ag films in 
the infrared is discussed in the literature [l|, |^, [ij, [H, Isil . 
m, [S^-Ull , summarized in Fig. [31 Selecting a refractive 
index out of these is an obvious source of large error in all 
results shown further down, and the largest uncertainty 
in modeling efforts of this kind j2| . Early design assumed 
values > 0.98 9, 30]. 

Because no ellipsometric information on any of our 
reflections is available, we start from a constant n = 
0.67 — 18. 3z for all surfaces, which is the Hagemann- 
Gudat-Kunz value of Silver interpolated to a wavelength 
of 2.48jim on logarithmic scales of the real and imaginary 
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part [22|. We ignore that Ml is clearly aluminized 
and that M9 is a dichroic. The net effect of this choice is 
an energy transmission coefficient (absolute value of the 
determinant of the Jones matrix or the Mjj element of 
the equivalent Mueller matrix) near 0.77. 

In a desperate act of tuning this to the anticipated 
transmission coefficient near 0.2 — geometrically scaled by 
mirror count to a different configuration 36] and remov- 
ing an estimated contribution from diffuse refiection — , 
this value of n is multiplied by 0.13. The calculations 
actually use n = 0.0871 — 2.379i. The average inten- 
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FIG. 3: Middle and top: Real part and absolute value of 
the i mag inar y part of measured refractive indices n of Silver 
0, S, [ia, El [2I, lli] as a function of wavelength A. Bot- 
tom: associated intensity reflection coefficient at zero inci- 
dence, \R\^ = \{n-l)/{n + l)f. 



sity reflection coefficient at the average incidence angle 
is brought down to \/0.2 0.95. 



B. Jones Calculus 

The ray tracing is started with a minimum distance to 
the optical axis, as enforced by the M2 diameter. A small 
corrugation of polarization across the beam as a function 
of the incidence point on Ml is induced by the variation 
of incidence angles on the mirror curvatures 0, [s^ E^l . 
It is not studied here. 

The Fresnel coefficients or each reflection are the two 
entries on the diagonal of each Jones matrix [l^ [2^, 



1271 . |35|, in the s— p— basis for the two states of polar- 
ization. The contribution of an individual reflection is 
the product with the 2x2 rotation matrix that adapts 
the plane of incidence of the previous mirror surface to 
the plane of incidence at the surface point of the next 
mirror, and with the free flight propagation term which 
is a diagonal matrix with two phasors depending on the 
(optical) path length of the ray from one mirror surface 
to the next [li,!!!]. 

The net effect of the mirror train train is the product 
of the matrices, four complex values J. It converts an 
initial Jones vector with amplitudes Ap_s into a vector 
with amplitudes Rp^s (Figure |4]), 
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FIG. 4: Illustration of the main ellipsometric variables that 
characterize the polarization state. 

The values plotted in figures [SH?] consider circularly 
polarized light with 



A., 



(2) 



as the input above Ml. The relative phase shift is 
S = Ss — Sp The ellipse main axis is rotated by an 
angle a relative to the main coordinates (where Ap is 
the horizontal component as we consider the beam in the 
interferometric laboratory where the planes of incidence 
are horizontal): 
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This formula defines < a < 90° as the direction of the 
longer or of the shorter axis. One can attach it uniquely 
to the longer axis if the sign of Rp is flipped and Sp shifted 
by 180° for all cases where cos 5 < 0. This will add 
90° steps in the foUowup graphs of a wherever the el- 
lipse passes through a "degenerate" circle as a function 
of pointing direction. 
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Let the amplitudes of the ehipse in its main coordinates 
be a and b (Fig. The amphtude (axis) ratio of the 
polarization ellipse then is 



tanx = -, 
a 

which is computed from a and S via @ 
sin(2x) = ± sin(2a) sinJ. 
The diattenuation 

62- 



P = 



&2 



= |cos(2x)| 



(4) 



(5) 



(6) 



measures the asymmetry of the ellipse, converted from 
field amplitudes to energy or flux. (In Figure SI b/a fa 
0.44 and therefore P « 0.67.) The value of P turns out 
to be the same for the left- and right-handed circulation 
of the circular polarization, that is, to be immune against 
flipping I to —I in the lower component in The value 
of a changes by 90° under this reversal of the input po- 
larization. 



III. RESULTS 
A. Diattenuation 

Figures [5] to [7] show the influence of the pointing di- 
rection (A, z) in the topocentric horizontal system and 
of the turning angle r of the D1-D3 triple on the axis 
ratio P and on the tilt angle a/° of initially circularly 
polarized star light after its turn into the +V direction 
inside the laboratory. Two plots of each figure deal with 
the N+ configuration, the two other plots with the N- 
configuration. 

In the plots of P{A,z), the pointing direction is en- 
coded in the circular coordinates of the base: the azimuth 
A = is indicated with a blank triangular sector. Point- 
ing to the zenith, z = 0-5°, happens in the middle of the 
small blank center. Five circular segments, each span- 
ning a zenith interval of 10°, lead to the outer rim of the 
plot where the zenith angle z reaches 55°. The value of P 
is plotted as a function over these {A, z) coordinates and 
in addition shown on color scales (dark blue for small P, 
light orange for large P) inside the base circle spanned 
by A and z. The value of a is shown as a scatter plot as 
a function of A for six different zenith angles z. 

Each of the figures starting with Fig. [5] remains un- 
changed if one switches between the -\-U and the —U 
position of the mirrors in the main delay line tunnel (mir- 
rors of the main delay line plus the associated M12 and 
M16 in the old nomenclature). 




N+, z=5 
N+, z=15 
N+, z=25 
N+, z=35 
N+, z=45 
N+, z=55 



, z=5 
z=15 
z=25 
z=35 
z=45 
z=55 



FIG. 5: The diattenuation ((6]) and angle a ([Sj in degrees as 
a function of telescope pointing at fixed rotator angle r = 0. 
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FIG. 6: The diattenuation ([6]) and angle a Q as a function of FIG. 7: The diattenuation (|6} and angle a Q as a function 
telescope pointing at fixed rotator angle r = 30° . At azimuths of telescope pointing at fixed rotator angle r = 60° . 
A for which P = 0, keeping a along the longer axis of the 
ellipse creates artificial jumps of 90° in a. 
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In overview, each of the two constehations N+ or N- 
of the Nasmyth station has two azimuth directions with 
preferred, low asymmetry P inside the beam, indicated 
by two blue-shaded radial sections on the base of the 3D 
plots. These directions co-rotate with the rotator angle 
r, which demonstrates that "in principle" this degree of 
operational freedom could be used to reduce polarization 
effects. In practice, the rotator angle r will often be de- 
ployed to serve other needs, which are not a topic of this 
manuscript. 

In Figure [5l the rotator mirror D2 is aligned with the 
other mirrors of the Coude train, and the main observa- 
tion is that the asymmetry P decreases while increasing 
the zenith angle — an effect of decreasing incidence angles 
on M3. 

The second, main observation is that P is large in the 
other two, unfavorable pointing directions A, whereas the 
dependence on the zenith angle is comparatively weak. 

Third, the difference between the A^-|- and N— config- 
urations for any prescribed pointing direction is smaller 
than one might have hoped, although there is some in- 
fluence on P. 



B. Mueller Matrix 

Another view on the same data is taken by transform- 
ing the Jones matrices into Mueller matrices [8., ,43] . 

Results referring to the Q and U components of the 
Stokes vector depend on the choice of the coordinate 
axes of the s- and p-polarization. Because the mirror 
train is axially symmetric with respect to reflections off 
Ml and M2, the p-direction above Ml has been defined 
relative to the plane of incidence on M3, the first mir- 
ror to break this symmetry. As we are co-rotating the 
Nasmyth focus with the azimuth pointing angle A, the 
Cartesian axes of the polarization directions above Ml 
are therefore defined not a in a coordinate system laid 
out by the architecture of the duct and main delay line, 
but in the coordinate system spanned by M2, M3 and 
M4. The direction of positive p is parallel to the direc- 
tion from MB to M4 of the N-|- configuration. The main 
intent of this comoving polarization frame in the input 
pupil is to avoid steps in the definition for azimuth an- 
gle transits through some fixed value; as a disadvantage, 
the representation contains an artificial smooth change 
of this frame through 360° for one full rotation around 
the vertical axis. 

The components in the interferometric laboratory re- 
fer to the natural horizontal coordinate system. (The 
unfortunate standard names U, V and W of the global 
coordinate system used in Figs. [TH2] have nothing to do 
with these assignments of polarization states.) 

Figure [5] indicates a high survival rate of star light 
polarization: the three relevant diagonal elements of the 
Mueller matrix have coupling amplitudes of the order of 
the intensity transmission of « 0.2. 

Figures [5HTil illustrate with the top row elements the 



conversion of the /, Q, U and V components of the 
Stokes Vector into an / component on exit of the mir- 
ror train. These quantities answer the question how far 
a preferred state of polarization of the beam above Ml 
induces changes in intensity observed in the interfero- 
metric laboratory. Their top graph is Af//, the intensity- 
to-intensity conversion by the mirror train, which was 
forced to a value near 0.2 by the scaling explained in 
Section Hi Al The values of the off-diagonal Mjq, Mjjj 
and Mjv in the other graphs vary typically over a range 
±0.01, which indicates that fully polarized star light may 
lead to relative intensity variations of up to 10%. 

As already apparent in Figs. [SHZl positioning of the 
rotator in its "neutral" direction r = leads to the least 
sensitivity of the matrix elements on the pointing direc- 
tion: The spread of the Af- values in Figs.[5HTUlis smaller 
than the spread in Figs. [TTHT41 




FIG. 8: Three diagonal elements of the Mueller matrix as a 
function of telescope pointing at fixed rotator angle r = 0°, 
Nasmyth configuration N-I-. 
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FIG. 9: The four elements of the top row of the Mueller FIG. 10: The elements of the top row of the Mueller matrix 
matrix as a function of telescope pointing at fixed rotator at fixed rotator angle r — 0° , Nasmyth configuration N-. 
angle r = 0°, Nasmyth configuration N+. 
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FIG. 11: The elements of the top row of the Mueller matrix FIG. 12: The elements of the top row of the Mueller matrix 
at fixed rotator angle r — 30°, Nasmyth configuration N+. at fixed rotator angle r — 30°, Nasmyth configuration N-. 
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IV. SUMMARY 

Apparently this is the first publication to assess the 
VLTI beam polarimetry induced by reflections off the 31 
mirrors of a standard optical configuration. 

The results remain semi-quantitative because they are 
based on a blind estimate of the dielectric function of 
the mirror surfaces, so only the dependence on the vari- 
able parts of the mirror train geometry has been empha- 
sized. To first order, interferometry with a symmetric 



setup, sensitive to differential properties of two telescope 
beams, is not affected, and the calibration procedure se- 
lecting calibrator sources nearby the science target will 
wipe out most of the pointing dependencies that were 
outlined here. 

For circularly polarized star light, the diattenuation 

(visibility contrast tested with a rotating analyzer) of the 
beam is predicted to reach values up to 1.0 for unlucky 
pointing azimuths, which are themselves a function of the 
star rotator position. 
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